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^ Abstract 

• ■ The Curie- Weiss Potts model is a mean field version of the well-known Potts model. In 

' this model, the critical line (3 = Pc{h) is explicitly known and corresponds to a first order 

, transition when q > 2. In the present paper we describe the fluctuations of the density vector 

in the whole domain /3 ^ and h ^ 0, including the conditional fluctuations on the critical 
line and the non-Gaussian fluctuations at the extremity of the critical line. The probabilities 
of each of the two thermodynamically stable states on the critical line are also computed. 
Similar results are inferred for the Random- Cluster model on the complete graph. 

> 

in 

^ ■ 1 Introduction 

The Curie- Weiss Potts model is a model of statistical mechanics which, being a mean-field model, 
can be studied by means of analytic tools. First it was shown in [1] that at h — Q, the model 
undergoes a phase transition at the critical inverse temperature 



OO 

O ■ _ / 9 if 9 ^ 2 

^^-1 2f3ilog(g-l) ifg>2. 



>< 



When q > 2 this transition is first order. The case of non-zero external field was considered in [2] 



^ I and it appeared that the first-order transition remains on a critical line. Recently this critical line 

was computed explicitly [3]. 

On the critical line, two or more states can coexist. One of the issue we address in the present 
work is the computation of the probabilities of these stable states. We also obtain a description 
of the limit distribution of the empirical vector of the spin variables that extend previous results 
on the Curie- Weiss Ising model [Ij (see also [HE]), and previous results on the Curie- Weiss Potts 
model with no external field [7]. 

The Curie- Weiss Potts model is connected as well to the random-cluster model. In that model, 
the first order phase transition for q > 2 was described in [S] and it appeared that at criticality, 
two possible structures of the random graph are possible. The probability for each structure was 
latter computed in . A consequence of our results we present a simple way of computing these 
probabilities when g > 2 is integer. 
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2 The Curie- Weiss Potts model 

The Curie- Weiss Potts model is a spin model on the complete graph. The probability of observing 
the configuration cr G {1, . . . , g}" at inverse temperature /3, in an exterior field H = h/ (3 equals 

^i^3,h.nW) ^ IT- exp I - ^ <^o-,,o-,- + /ty^(5cr.,i 



I3,h,n 



where 5 is the Kronecker symbol and Z^.^.n the partition function 



■'PJi.n 




Our interest is in the limit distribution of the empirical vector 

Cn n \ 

(2-1) 
1 = 1 4=1 / 

that represents the number of spins of each color for a given configuration a. The normalized 
vector N/n belongs to the set of probability vectors 

{x eW -.xi^ h a;, = 1 and ^ 0,Vi}. (2.2) 

The large deviation principle for N/n is an immediate application of Stirling's formula (see for 
instance Lemma l4.ip . If we consider f^^h the microcanonical free energy of the model: 

fpji{x) ^^XilogXi - -^Xi - hxi, Vxen (2.3) 

i=l 1=1 

with the convention that log = 0, then we have the following classical large deviation result (see 
for instance [TD], and also [SHU] for LDP concerning the closely related random cluster model). 

Theorem 2.1 Assume that /3„ — > /3 and hn — *■ h. Then, the vector N/n £ 51 distributed according 
to the measure fJ,f3„^h,^,n follows a large deviation principle with speed n and good rate function 
fp^h - mino fj},h. 

This large deviation principle leads to a law of large number: when ffj^h has a unique global 
minimizer, N/n converges towards that minimizer. The structure of the minimizers of ffj^h was 
determined in the papers [ll[l2j[2j. Here we give some further details: 

Proposition 2.2 Let P,h ^ and let x be a global minimizer of fp_h in ^■ 
i. The vector x has the coordinate min(2:i) repeated q ^ I times at least, 
a. If h> Q, then xi > Xi, for all i G {2, . . .,q}. 
Hi. The inequality min(a;i) > holds. 

iv. For any q ^ 3, or q — 2 and (/3, h) ^ [(3c, 0), one has min(xi) < 1//3. 

Because of the simple structure of the global minimizers of the free energy, the problem of finding 
them reduces to a one-dimensional optimization problem. The usual parametrization consists in 
taking xi — [1 + {q — l)s)/q, X2 = ■ ■ ■ = Xq = [1 — s) / q where s £ [0, 1] is a parameter called the 
magnetization. Another equivalent parametrization permitted in [3] the explicit computation of 
the critical line 

/it = |(/3, /i) : ^ /i< /lo and = log(9 - 1) - Z^^^^j (2.4) 
with extremities (/?c,0) and (/3o,/io), where 

/3o = 4^— i and ha = log(q - 1) - 2^-^ 
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were already determined in [5]. The key observation in 3J was that the free energy fp^hixz) at 

'1 + z 1- z 1- z ' 



2 2(9-1)' 
is easily split into its even and odd parts: 



2(g-l) 



z e [±l] 



(2.5) 



1 + 1 + z I - z 

■ log — ! 7, — log ■ 



2 "2 

log(g - 1) - /5 



2 ° 2 
9-2 , 



ilog(,-l)-^^ii±^ 
2 Bvy y 8 
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2(9-1) 



showing that, on the critical line hr, the free energy fp.hixz) is an even function of z. It is strictly 
convex for f3 < (3o but not for f3 ^ (3o- Indeed, the second derivative of z i-^ f/3.h{xz) is 



dz2 



1 



[iq 



l-z2 4(9-1)' 



(2.6) 



thus, for /3 ^ /3o, the function is strictly convex on [— 1,— z^) and on (2:^,1], concave on {—Zi,Zi) 

where 

z, = v/1 - Pol 13- (2.7) 

Depending on the parameters (/3, h) the free energy presents one or several global mininiizers. The 
following is a summary of the works [T] (for h = 0) and [3^ (for /i > 0): 

Theorem 2.3 Let f3,h^0. 

i. If h > and {/3,h) ^ hr, the free energy fp_h has a unique global minimizer in fi. This 
minimizer is analytic in f3 and h outside of Ht U {(/3o,/io)}- 

ii. If h > and {f3, h) g Ht, the free energy fp^h has two global minimizers in O. More precisely, 
for any z G (0, {q — 2)/g), the two global minimizers of fp_,^h^ at 

Pz = 2^—^- log and = \og{q - 1) - 77^^^ — '^Pz 



q z 



2(9-1) 



are the points x±z- Furthermore, (resp. x^z) is the limit of the unique global minimizer 
of ffi.h o-s {P,h) —f (Pzjhz) above (resp. below) the line hx- 

Hi. If h = and /3 < Pc, the unique global minimizer of fp_h is {1/q, . . . , l/f?) = a;_(g_2)/g- 

iv. If h ~ and P > Pc, there are q global minimizers of fp,h, which all equal Xz up to a 
permutation of the coordinates, for some appropriate z G ((g — 2)/(7, 1). 

V. Ifh = and P — Pc, there are q+l global minimizers of ffj^h ■' the symmetric one (1/q, . . . , 1/q) = 
a;_(g_2)/(} together with the permutations of 



9-1 1 
9 '9(9-1)' 



9(9-1; 



= X 



(9-2)/q- 



3 Statement of the results 

In this paper we address essentially two questions. According to Theorem 12.11 the distribution of 
N /n is concentrated, as n ^ +00, on the set of global minimizers of the free energy. First, we 
study the fluctuations of the empirical vector TV around its typical value. Second, when several 
global minimizers exists we explicit the weight of each of them. 

These questions were answered in several very interesting papers for particular cases of the 
model. The case of the Curie- Weiss Ising model {q — 2) was reported in 0] (see O [6] for the 
proofs), while the Curie- Weiss Potts model was treated at zero external field in [7]. 

Our approach is similar to that of the former references, with the technical difference that our 
computations are based on Stirling's formula while the former works are based on the fact that the 
law of N/n + W / y/n, where is a Gaussian vector in R'' with distribution Af{0, P~'^Iq), can be 
explicitly computed (see for instance Lemma 3.2 in [7]). 
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We also permit that the parameters [3 and h fluctuate with n, and take in the sequel /i„) — > 
{P,h). This will be useful for applying our results to related model such as the random cluster 
model on the complete graph. 

Our first result concerns the fluctuations of the empirical vector TV outside of the critical line. 
The fluctuations belong to the hyperplane 



n = jtt; e R'' : = o| . (3.1) 

Not surprisingly, these fluctuations are Gaussian. This generalizes Theorem 2.4 in to the case 
of positive external fields. The way that (/3„, hn) converges to {/3, h) is able to shift the center of 
the distribution. 

Theorem 3.1 Assume that (/9„,/i„) {I3,h) for some j3,h^Q with {13, h) ^ (/3o,^o)- Assume 
that there is a unique global minimizer x = (xi, Xq, . . . , Xq) of the free energy fp,h- For every n, let 
dn the smallest d G Ti such that x + d Cz is a local minimizer of fp^^h^. Let W be the random 
variable in Ti such that 

N = jix + nd,-, + n^/'^W (3.2) 

where the distribution of N is given according to the measure fJ-f3„^hn,n- Then, W converges in law 
towards the centered Gaussian vector with covariance matrix 

( q-1 -1 ••• -1 \ 



q(3 



-1 l + iq-2)- 



l + (.-2)|^y 



(3.3) 



which has rank q — 1. 



Remark 3.2 The vector <i„ is O (|/?„ — (3\ + \hn — h\) when the quadratic term in the Taylor ex- 
pansion of fp^h is definite, that is for any (/?, h) ^ (/3o, ho) - see Lemmas \4.5\ |^.6] and \4.7\ below. 
Hence, for Pn — f3 — o{n^^^^) and hn ~ h — o{n~^^'^) the vector ndn is negligible with respect to 
n^f^W and could be removed from the definition of W at L3.2\) . It is remarkable also that on the 
line (3 < [3c, h — we have x ~ {1/q, . . . ,1/q), hence for /i„ ~ and /3„ (3, the vector <i„ is 
exactly zero. 

Remark 3.3 In the range of validity of Theorem 2.4 in ^Tj, that is (3n = P < Pc cind hn = h = 0, 
we have xi — . . . — Xq — 1/q thus the covariance matrix simplifies to 



1 



/ q-1 -1 



q^ - qP 



V -1 q-l 

We have checked the correspondence with the covariance matrix that appears in J^. 

The matrix (|3.3p gives a special emphasis on the first coordinate since it corresponds to the 
case X2 = . ■ . = Xq. Before stating the next theorem we give a more symmetric definition for the 
covariance matrix: we let 




-a{xi,Xq) 

1 + (g - 2)a{Xq,Xq) 



\mm.[Xi)ma^(Xi) J 

where 

a{x,y) 



max(a;. y) ^ ~ P 

When the free energy has several global minimizers, that is when (/?, h) G hx or P ^ Pc and 
h = 0, the empirical vector 7V/n is close to either one or the other of the minimizers of the free 
energy ff^^h- We first determine the conditional fluctuations (this extends Theorem 2.5 of [7]): 
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Theorem 3.4 Assume that (/?„, /i„) — > h) with f3,h ^ 0. Assume that the free energy fp^h has 
multiple global minimizers x, x' , . . . and let s > smaller than the distance between any two global 
m,inimizers of fp^h- Let dn the smallest d G 7i such that x + d € ^ is a local minimizer of fp^.h^ ■ 
Then, under the conditional measure 



— e B(x,e) 
n 



the variable W defined by N — nx + nd„ + n^/'^W converges in law to the centered Gaussian 
vector with covariance matrix K{x). 

Additionally we compute the limit probabilities that N /nhe close to a given global minimizer 
of the free energy, generalizing Theorem 2.3 of fl\ with an explicit formula. 

Theorem 3.5 Assume that there are P,h ^ and A, G IR such that 

(/?„, h„) - h) + (A, ly) + o{n-'), 

and assume that the free energy fp^h has multiple global minimizers x,x'^... If e > is smaller 
than the distance between any two global minimizers of fp,h, then 

/N \ t(x) 

lim /./3„,„„,„ — GB(a;,£) = ^ ^ ^ (3-5) 



wht 




t{x) = ( 1 — /? min(xi) ) exp I ' xf + vxi I . (3.6) 



Remark 3.6 On the critical line hx one can parametrize the formula i fg.5)) according to the second 
point in Theorem \ 2.3\ : when {f3,h) = (/3z,/iz) with z G (0, (g — 2)/q) the two global minimizers 
are x±z. In particular, when z — > (i.e. {P,h) on hx close to (/3o,/io)J, the probability of each 
corresponding state converges to 1/2. 

We also describe the fluctuations at the extremity (/3o, hg) of the critical line. This extends for 
instance Theorem 2 in [4] that applies to the case of the Curie- Weiss Ising model at criticality, 
namely q = 2 and (/3o,/io) = (/3c, 0). We recall that H defined at p.ip is the hyperplane parallel 
to f2. Given a vector m G R"^, we denote by the vector space made of all vectors orthogonal to 
u in the Euclidean space R'. 

Theorem 3.7 Assume that (/3„, ft,„) — > (/3o, ha) with /3„ — /3 — o{n^^^'^) and hn — h = o{n^^/^) and 
let X = (1/2, l/2((7 — 1), . . . , l/2((j — 1)) be the unique minimizer of fp^^ho ■ Let m = (1 — g, 1,...,1). 
// the random variables T G IR and V G 7i H u-^ are defined by 

N = nx + n^l'^Tu + n^l^V, (3.7) 

then (r, V) converges in law. The limit has the following properties: 

i. T and V are asymptotically independent 

a. T converges in law to the probability measure on IR proportional to 

e.p(-M,.)d, 

Hi. V converges in law towards the centered Gaussian vector with covariance matrix 

/ ••• \ 

q-2 -1 



2(g-l)2(g-2) 

of rank q — 2. 



\0 -1 q-2 J 
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We conclude the summary of our results with two claims on the random-cluster model G{n,p, q) 
on the complete graph Kn with n vertices. In that model, a configuration lj e {0, has a 
probability proportional to 

where C{lo) stands for the number of connected components of the sub-graph with edge set {e G 
E{Kn) : We = 1}. This model is closely related to the Potts model after the well known Fortuin- 
Kasteleyn representation (see for instance [13]). We take a spin configuration a G {1, . . . , 
under the measure fJ.[3„^h„,n, then let = 1 with probability pn — 1 — exp(— /3„/n) only if ct; — Uj, 
where i,j are the extremities of the edge e (and else Ue = 0). The resulting configuration uj follows 
the distribution of the random cluster model G{n,pn, q)- 

First we have a Corollary of Theorem 13. 51 we compute the probability that there exists a giant 
component in G{n,Pmq), that is a connected component for u of size Q{n), when p„ is close to 
the critical value (5c/ n. This completes part (b) of Theorem 2.3 in [5], with a simpler proof than 
that of Theorem 19 of 0. 

Corollary 3.8 Let q> 2 integer and consider pn such that 

Pc^l ^ (I 
Pn^ ^ ^+°[^ 

Then, with a probability that converges to 

1 



i+KT-i^)"»-(-(?+i)Sf^ 

the graph G{n,pn, q) contains a giant component. 

The description of the Gaussian fluctuations also enable fine computations of the partition 
function of the random-cluster model 

Z^ln = E n P"=(1-P)^-"='Z'"("^. (3.8) 

For instance, 

Proposition 3.9 The partition function of the random cluster model for integer q ^ 2 and 

/3 , 7 , /I 



Pn — I 9 ' " I 9 



with ^ P < Pc ttnd 7 G R satisfies 



P\ ^ n f 2nf3 + 2j + pWq-l 



Z^^,,,n = {l + o„{l))[l-^) g"exp(^- "-" ' 2' [^))- ^^-^^ 

Remark 3.10 Although our Theorem \3.5\ agrees with Theorem 2.3 of ^ when h — 0, Corollary 
3.8\ and Provosition \3. 9\ do not give exactly the same conclusions as, respectively. Theorem 19 and 
Theorem 9 (i) in l^j. The latter Theorem states an equivalent to the partition function restricted to 
the set of configurations made of trees and unicyclic components, which, for P < (3c, is equivalent 
to the whole partition function. The ratio of the equivalent in Theorem 9 (i) in JS/ over ^3. 9\) is 



(the formulas do coincide at the exponential order). We could not check the proofs in 191, yet we 



were surprised to find that Theorem 9 (i) would not permit to recover Z^^i „ = 1 for q ^ 1. 
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4 Proofs 



This section is organized as follows. First we describe the asymptotics of the distribution using 
Stirling's formula. We also prove Proposition l2.2l Then we address successively the limit distribu- 
tion at h) ^ (/So, ho), at the extremity (/3o, ho) of the critical line, and finally we give the proofs 
of the related results for the random-cluster model. 

4.1 Asymptotic density & limit of the uniform measure 

In this Section we give an equivalent to the density of the Potts model, prove Proposition 12.21 and 
describe the limit of the uniform measure on the set of possible realizations of N/n. 
For any e ^ 0, we let 

fig = {a? e : minxi ^ e} and il" = {x e fls ■ nx e N'}, 

where f2 is the set of probability vectors, see ()2.2p . We also write ilQ+ = {J^^Q^te and fl" = ^Iq. 
In our first Lemma we give an equivalent to the density of the Potts model with respect to the 
counting measure on fi". We use nothing else than Stirling's formula 

nl = (l + o„(l))V2^(-j . 

For any x S r2o+ ^-nd /3 ^ 0, we let 

Af,ix) = (27r)-^ n^r'^'exp . (4.1) 

We also recall that the free energy fi3,h was defined at (|2.3p . We have: 
Lemma 4.1 For any f3,h, any 1 and x £ define rfj,h,n{x) by 

Zi3,h.nfJ-i3.h,n{N ^ nx) = {1 + r p ,h^n{x))n ~ Ap{x) exp{-nfpji{x)). 
Then, for any e > 0, supj.gQ„ sup^_/j \'r'/3^h,n{x)\ goes to as n ^ oo. 

Proof Given x e 17" we write nx — (ni, . . . , ng) n. It is a vector with positive integer 
coordinates. There are exactly nl/ Y[i=i ""^i ways of choosing the spin configuration that satisfy the 
constraint (A'^i, . . . , Nq) — (ni, . . . , n^), hence 

, , n\ IP ni(ni — 1) \ 



i=l 



ni! ■ ■ -Uql 



■ exp \ —— + n 



2 

i=l 



Thus 

1 + rp^h.nix) 



Zl3,h,nH0,h,n {N = nx) 

n~^Ap{x) exp{-nfpji{x)) 



nil ■ ■ ■ Uql 



. exp [n} log Xi 

V27rn \ J 



which does not depend on /? nor on h. Applying Stirling's formula yields the conclusion as all the 
Ui go to infinity uniformly over a; e fi". □ 

Remark 4.2 Theorem \2.1\ is a consequence of Lemma \4-l\ as fp,^,h„{x) — * ffs,hix) when (/3„, ft,„) - 
{f3,h), uniformly over x G (for the uniformity, see Lemma \4-.6\ below). 
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Now we give the proof of Proposition 12.21 For completeness we repeat some arguments from 

[mi. 

Proof f Proposition 12 . 2p . Call gp{z) — log(z) — i3z. As a; is a minimizer of the free energy in H 



one has 



oxi oxj 



gpixt) - 9i3{xj) - h[li,=iy -ly^iy] , yi,j G {l,...q} (4.2) 
and < ^(-) + ^(-) 

9'f3{xi)+9'p{xj)- (4.3) 



First we assmne h — 0. As 5 is concave, ()4.2p implies that the set {xi : i — 1, . . . ,q} contains at 
most two values. Equation (|4.3|) imphes that at most one of the Xi has g'p{xj) < 0. As 

' X 

is positive on (0, 1//3) and negative on (1//3, 1), the first point of Proposition 12.21 follows together 
with the inequality 

min(x,) s: 4- (4-4) 

P 

Assume now that ft, > and that (|4.4p does not hold for some i E {2, . . . ,q}. If Xi > xi with 
i G {2,..., 5}, the vector x with xi and Xi permuted has fp^h{x) < ffjji{x), a contradiction, 
therefore xi ^ a;,;. The equality xi — X2 is impossible in view of (|4.2p . yielding the second point of 
Proposition Now wc conclude the proof of the first point of Proposition 12 . 21 when h > : the 
inequality Xi ^ 1//?, that implies g'pixi) ^ and g'pixi) < since xi > Xi, would contradict (|4.3p . 
Hence all the Xi belong to (0, 1//3) where there is at most one reciprocal image of gfsixi) — ft, by 
gis, hence X2 = ■ ■ ■ = Xq < 1/ f3. 

Now we address the third point. If min(a;i) = 0, one can find i,j £ {1, ...,(?} such that Xi = 
and Xj > a.s X £ 57. Hence x* — x + t{ei — ej) belongs to ilo+ for small enough t > 0. Yet, 

^//3,/i(a;*) = gpixi+t)- gpix-i -t) + h[l{i^iy-l[j^iy] 

goes to— ooast-^0+, a contradiction. 

Remains the strict inequality in (|4.4p . We let x^^ = (1 — {q — l)/f3, . . . , 1//3) (which is in 
for /? ^ g — 1, and satisfies the case of equality in (14. 4p for f3 ^ q) and derive conditions for X/^ 
being a minimizer of the free energy. Equation (j4.2p for i = I, j = 2 gives 



ft = log(/3a;i) — Pxi + 1 

which is negative unless xi also equals 1/(3. Yet, xi = implies (3 — q and ft = 0. But 
a; = (1/(7, . . . , 1/g) is a minimizer of the free energy fp^ only for /3 ^ /3c- As g > 2 ^ /3c < the 
only case of equality is g = /? = 2 and ft = 0. □ 

In a second Lemma we compare the counting measure on 57" with the Lebesgue measure. This 
will help in the proofs of Theorems 13.11 13.41 13.51 and 13.71 We denote by £ the Lebesgue measure 
on hyperplanes. 

Lemma 4.3 Let H : E,'' ^ R"* be a affine and one-to-one transformation. Let P = [0, 1]'' n H.. 
Then, for any f : IR' — s- IR bounded, 

where 

ip{z) = sup f, VzeR«. 

n-i(z)-n-ip 
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Remark 4.4 Applying this to — / one obtains a useful lower bound. 

Proof For any 2; G 11 (a; + n^^P) , one has 11^^(2;) e a; + n^^P hence x G 11^^(2) — n^^P. Thus 
tf{z) ^ /(a;), and 

C{\l{n ip)) 7n(x+n-ip) 

The claim foUows when we sum over x G as il" + n^^P = f2 + n^^P. □ 
4.2 Gaussian fluctuations 

The hmit theorems will be proved as consequences of a Taylor expansion of the free energy. First we 
consider a second order expansion of //3,/i, that will be enough to describe the Gaussian fluctuations 
at (/3,/i) ^ (/?o,/io). 

This section is organized as follows. First we give a series of Lemmas that permit to establish 
Proposition [TS] below. Then we give the proofs of Theorems 13.11 and 

4.2.1 Taylor expansion of the free energy 

The Taylor-Lagrange formula applied to the C°° function t G [0, 1] 1-^ fp,h{x + tw) yields: 

Lemma 4.5 Let x (z ^ be a global minimizer of fp^h o,nd w ^TL such that a; + it) G Oq+ . Then, 
there exists a G (0, 1) such that 

fpAx + w) = /^^^(a;) + i^f— ^ Awl (4.5) 

On the other hand, the influence of /?„ and /i„ is immediate to characterize: 
Lemma 4.6 For any f3, /3„, h, h„ and any a; G il, the following equality holds: 

fp^.hAx) = fl3,h{x) - ^"^ ^ ^xj - [hn - K)xi. (4.6) 

1=1 

4.2.2 The quadratic form 

Given x G Vta+ and /3 ^ we consider the quadratic form Qx.fi : s- IR defined by 

9 



i—1 ^ * / 



(4.7) 



This is the quadratic form that appears in Lemma 14.51 When it is positive definite it determines 
the fluctuations. We have: 

Lemma 4.7 Let f3,h ^ and x be some global minimizer of fp^h- 
i. The quadratic form Qx,p is positive definite on TL if and only (/3, h) 7^ (/3o, ft-o)- 
ii. When {f3,h) — (/3o,/io) o.i^'d x = Xq, the kernel ofQx^p is Vect(it) where 

u = (1 - q,l,. . . ,1) ■ 

Proof First we assume that (/3, h) ^ (/Jq, Hq) and prove that Qx,f3 is positive definite. According 
to Proposition 12 .21 the vector x has one coordinate repeated at least q—1 times. Let j G {1, . . . ,q} 
be the smallest index such that Xj — max(a;i), and J ^ {1, . . . ,q} \ {j}. For any w E H one has 
Wj = - J2ieJ hence 



-"x(a;,) J J \mm{x^) J ^ 
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Now we let 

, 2 



9-1 E^eJ'^^ 

which belongs to the interval [0, 1] according to Cauchy-Schwarz inequality, and obtain 



Qx{w) = 



max(xi) J \min(xi) 



J2wl (4.8) 



ie.J 



Hence the quadratic form is positive definite on Ti if and only if the factor in (|4.8p is strictly 
positive at both a = and a = I, that is to say if 

1 



min(a;i) 



13 > (4.9) 



and ^ , \ + .\ > 0. (4.10) 

max(Xi) mm[Xi) 

Condition ()4.9p is true as {/3,h) ^ {(3o,ho), cf. Proposition [^1 Condition ()4.10p is equivalent to 

1 — q(3mm{xi)m8iK{xi) > 

as (g — 1) min(xi) + max(a;i) = 1. The reader will remark that min(a;i) max(a;i) is constant over all 
the minimizers of fp^h described in Theorem 12.31 Hence we might take the minimizer of the form 
X ^ Xz as in (|2.5p . that is xi = (1 + z)/2 and X2 ^ ■ ■ ■ = Xq = {1 — z) / {q — 1) , which reveals that 
condition (|4.10p is equivalent to z i-^ ffi.h{xz) having a positive second derivative at its minima, 
as 

<^^ff3,h{x,.) ^ 1 (3q 

dz2 l_z2 4(g-l)' 

which is the case again as (/3, h) ^ (/3o, /io) (see the discussion after (|2.6p ). 

Assume now that {(3,h) = {(3o,ho). If q = 2, the quadratic form Qx,[3 is identically zero on 
H = Vect(u). li q ^ 3, we have h = ho > hence j — 1. The quadratic form vanishes at 
w ^ H ii and only if ai{w) = 1. In view of the definition of a, this is the case of equality in the 
Cauchy-Schwarz inequality: ai{w) = l4^W2 = -- -=Wq4^we Vect(M). □ 

4.2.3 Centering of the fluctuations 

As in Theorems 13 . 1 1 and 13 . 41 we let <i„ the smallest d G H such that x + d E £7 is a global minimizer 
of fi3„,h„- We have: 

Proposition 4.8 Assume that {Pn, hn) (/3, h) and let x E be a global minimizer of fp^h- 
i. For any R > 0, 

nfp^M^lx + dn + n'^/^w'j = nfp^,h„{x + dn) + ^Q^p{w)+On{l) (4.11) 
uniformly over w G Ti. H B{0, R). 

ii. If {(3,h) ^ {Po,h()), for small enough e > there is X > such that, for n large enough and 
any w £H with \\w\\ ^ en^^^, 

nfp„,h^(^x + dn+n-^/^w^ ^ nfp^j,^{x + dn) + X\\w\\^. (4.12) 
Proof We begin with an application of Lemma l4.5l at the global minimum point x + dn'. 

nfp^,h„ (x + dn + n-^/^w^ = nfp^j,^ [x + d„) + ^Q^+d„+a„-i/2^,/3(«^) 

for some a g (0, 1) depending on n and w. For ()4.1ip we only have to notice that d„ + an^^^'^w = 
o(l). For (14.121) we remark that as d„ ~> (cf. Remark [3.2p . for all n large enough and ^ en^/^. 



dn + an ^^'^w 



^ 2e. 



As this can be made arbitrary small, for small enough e the quadratic form Qx+d„+an-^/^w,p 
dominates ^/2, which is definite positive after Lemma l4.7l □ 
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4.2.4 Some linear algebra 



The next Lemma will be useful at the time of computing inverses or determinants. Denote by /„ 
the n X 71 unitary matrix and by An the n x n matrix with all entries equal to 1 . 

Lemma 4.9 Let M = aAn + hin- 

i. The determinant of M is 

det (M) = (6 + na) . 

(a + b)-^ ifn = l 

Proof We prove the first point as follows: let P(A) = det(A/„ — An) be the characteristic 
polynomial for the matrix An- The matrix An has rank 1 and eigenvalues 0, . . . ,0, n, which are 
the roots of the unitary polynomial P, thus P{X) = A"~^(A — n). The second point follows from 
an immediate computation. □ 



ii. When M is invertihle, it has 



4.2.5 Proof of Theorems [HH [an and 

As a consequence of Proposition 14.81 we give the proof of Theorems 13.11 13.41 and 13.51 

Proof (Theorem I3.ip . First we condition W = n~^/^{N — nx — ndn) on ||W^|| < R for some 

positive R. For g : IR' i— > R continuous bounded, Lemma [4. II and (|4.1ip in Proposition 14.81 vield 

JV/neSl" 

The transformation X = N /n ^ n(X) = W is affine. The image of by IT is greater than 
'Hr)B{Q, R) for large enough n as a; e 

rJo+, and on the other hand n(X + n-ip) = n{X) + n-^^^P, 
that is to say the dimensions of the image of the lattice element P go to zero. Hence Lemma 
gives 

Zp^.h„,nti0^.h„.n{g{W)l{^^w\Kn.}) = (l + o„(l))A^(a;)e-"/"-'-'^(^+'^") x 



1 



CiP) 



5(t(;)e"^^-.'^(™M£(i(;). (4.13) 

nnB{o,R} 



In other words the law of W conditioned on || ^ R converges to the distribution on 7Yni3(0, R) 
with density proportional to 

e"^'3-.'*(™) (4.14) 

with respect to the Lebesgue measure onHC] B{0,R). 

Now we show that the variable W is tight. Thanks to Theorem 12.11 we know already that for 
any e > 0, 

\hnsup ^f3^^^hn,n {\\N ~ nx\\ ^ en) — 0. 

n 

Thus it is enough to show that, for small enough e > 0, 

lini limsup/X0„,,,„,„ (||VF|| ^ sC en^/^) = 0. (4.15) 



According to Lemma |4T] and Proposition [48l for small enough e > and large enough k there is 
A > such that 



MA 



,h^,n{\\W\\^K\\\W\\^en'/^) ^ 



/^/3„,/>„,„ \\W\\ ^£^1/2) 



^ (l + o„(l)) 



2^Ar/nGO":K^||W|Keni/2 6 



E 



Ar/neO":||W||<, 



. e 



LvRm ^e--^ll"'ll'd/:(u;) 
^ (l + o„(l))- ^ ' 
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after Lemma Since Qx,i3 is positive definite the ratio goes to as k — > oo, giving (I4.15p . 

Let us show that this hmit distribution is as well the distribution of the centered Gaussian 
vector with covariance matrix (|3.3p . We take V a random vector in Ti. with the density (|4.14p with 
respect to the Lebesgue measure on Ti.. The law of V is also proportional to 

g-^Q..3('')d„2...dw„. 



This density can be expressed only in terms of the truncated vector V — (V2, . 
we take 



H 



1 



Xi 



1 



,Vq). Indeed, if 
(4.16) 



we have Qx,p{v) = *vHv and thus the covariance matrix of V is 



H- 



X /, 



9-1 




according to Lemma 14.91 Using the relation Vi — ~ X)i=2 ^ compute the remaining covariance 
coefficients, leading to the completed matrix p.3p . The rank of the matrix is not less that that of 
H, that is (J — 1, and it is also strictly less than q because of the linear constraint V E H. □ 

Proof f Theorem 13. 4p . The former proof can be repeated almost verbatim. One has to take care 
however that x needs not be of the particular form X2 — ■ ■ ■ ^ Xq (although it still has a coordinate 
repeated q~l times), and that the variable which is tight is W conditioned on N/n £ B{x, e). □ 

Proof (Theorem 13. 5p . The tightness of W conditioned on N/n S B{x,£) and the convergence 
of the law of W on bounded sets (cf. the proof of Theorem 13. ip imply that for any e > smaller 
than the distance between any two minimizcrs of //3,h, 



i^&B{x,e)) 



(l + o„(l)) 



A0(a;)e-"^''".''"(^+'*"(^»/, 



e 2 



dC{w) 



Hence we call 

C/3„,/»„,„(a;) = A^(a;)e-"/f'"'''"(=^+'*"(=^» / e-5Q.,3Md/:(u;) 



n 



and give an equivalent to Cf3^,h„.nix). First we compute the integral up to a constant factor. We 
pick j £ {1, . . . ,q} such that xj = max(a;i) and let J = {1, . . . ,q} \ {j}. The Lebesgue measure on 
Ti. is proportional to the measure induced on w by Jliej difi, given Wj — — X^ieJ (|4.16p 
we let 



max(a;i) 



P Aq-l 



min(a::i) 



P 1,-1 



and w = {wi)i^j, thus ^wHx.pw = Qx,i3{w) and therefore 

J "hi „ 7 



■q-l 



1 



min(a;i) 



q-2 



1 



min(xi) max(a;i) 



qP 



-1-1/2 



according to Lemma 14.91 If we multiply with the prefactor Ap{x) we obtain 



Ap{x) 



H 



'-{I- PwLUv{xi))^ (1 - g/?max(xi)min(x,)) ^^^4-17) 



as max(a;i) + (q — 1) min(a;i) = 1. Then we use Lemma 



"//3„,/i„ {X + dn{x)) 



^Pn__P_ ^ .(^3.-)-)2 _ ^^y^^ _ ^-j 1^,^^ ^ di{x)) 

i=l 



= nfp^h{x) - -^x1 -VXi+ On{l) 
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as dn{x) — 0{\/n) and a; is a global minimizer of ffs^h- Thus we have shown that 

dC f A /? 

C/3„,/i„,„(a;) = (l + o„(l))|= —expi-nf,3,h{x) + -'^x'f + i^xi-- 

2 — g 1/2 

x(l — /3niin(xi))~2~ (l _ g/3max(a;i) min(a;i)) 

where the factor dC/ Jliej does not depend on J. The claim follows from the remark that the 
product min(a;i) max(a;i) is constant over all the global minimizers x, x' , ... of the free energy /^^/j 
at any (/3, /i) e hr, cf. TheoremO □ 



4.3 Limit theorems at criticality 

The proof of Theorem 13.71 relies again on a Taylor expansion of the free energy: 

Lemma 4.10 Let h) = {/3q, Hq) and x = Xq = (1/2, l/2(q — 1), ... , l/2{q — 1)) be the unique 
minimizer of fp^h- Let m = (1— g, 1,...,1). For alltEM, and v G TiDu^ such that x+tu+v E ilQ+ , 
there are a, a' G (0, 1) such that 

4 1 4 

Uj 



1 t 

fpji{x + tU + v) = fl3,h{x) + -Qx+fu+av,l3{v) + 



Furthermore, 

(.^ (g-l)(g-2) „^,„ 

T^Q^.pH = \\v\\ 

z q 



i—l ^ 

Proof A second-order Taylor expansion in v yields 



fl3,h{x + tU + v) = f0Ji{x + tu)+\/fpj2{x + tu) ■V+^Q^+tu+av,p{v) 

for some a G (0, 1). The last q — 1 coordinates of the gradient Vfi3^h{x + tu) are equal, hence it is 
orthogonal to v. Then a fourth order expansion in t gives 



for some a' £ (0, 1). Indeed, the first order term is zero as x is the global minimizer of fp^h- The 
second order term is Qx,p {tu) /2 = in view of Lemma 14.71 Hence the third order term is 0, 
yielding (gHH]). 

Let us prove the last two formulas. The assumption v E HD implies vi = 0, hence 



1=2 



iE(2(.-l)-4i-l).? 

(g-l)(g-2) ^ 2 

q ^ 

^ 1=1 



On the other hand: 



g 4 



E ^ = 8(g - 1)4 + {q- l)8(g - 1)^ = 16(g - 1)4 



■ 1 



□ 



Using Lemma [4. 101 we establish the analog of Proposition [4 
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Proposition 4.11 Assume that (/3„,/i„) — > (/3o,/io) with Pn — Po — o{n ^Z'*) and /i„ — ho — 
o(n~^^'^), and let Xq £ fl be the unique global minimizer of fp^^ho- 

i. For any R > 0, 

ii^Ilii^||.|P + l(,-l)-V + o„(l) (4.19) 
q 3 

uniformly over v E HO u-^ O B{Q, R) and t E [—R, R] ■ 

a. For small enough e > and large enough R, for n large enough, for any v E Ti, O , 
i e ]R \ [-i?, R] such that \\n-^^Hu + n-'^/'^v\\ ^ e, 

nff3„ji„ (x + n^^/'^tu + nr^/'^v^ > nfp^h (x) - n ^" ^ ^° ||a;|p - n(/t„ - ho)xi 

+ \\v\\yA + t^/2 (4.20) 

Proof We first apply Lemma for w = n^^^^tu + n^^^^v, we have 

"■//3„./i„ (a? + i«) = nfi3^h{x + w)-n-^ — ^^||cc|p - n(/i„ - /io)a;i + o„(l) 

uniformly over v E Hr\u^r\B{Q, R) and t E [-R, R] as P„-f3o = o{n^^/'^) and h„-ho^ o{n^^/^). 
Then, Lemma [4.101 vields 

nfp^,,{x + w) = nffj,h{x) + ^^^^^^^^\\vf + Uq-^ft' + On{l) 

q o 

uniformly over the same domain, and (j4.19p follows. 

Now we only assume that \\w\\ ^ e. For small enough e > we have, after Lemma f4.101 the 
lower bound 

nff}.h{x + w) ^ nff3,h{x) + ^\\vf + t^ (4.21) 

(note that, for g = 2, is necessarily 0). Combining with Lemma 14.61 we obtain that, whenever 
Km holds, 

nfi3^,hAx + w) ^ cnix) + \\v\\^ ( ^ ~ ^" ~ ^° ) - Int - m^t^ + 



where 



Cn{x) = nfi3ji{x) - n ^" ^ ||a;|p - n(/i,„ - ho)xi, 
In = n''/Hfin-Po)x-u-n^'\hn'ho){q-l) 

l/2/3n — /3o II II 2 

m„ = n ' ||m|| 

a.s x,u J- V and vi = 0. Now we conclude: for any large n, 



Similarly, as Z„ = o„(l) and m„ = o(n ^/^) ~ o„(l), for any large n and t large, 

t^ 



Int - rUnt^ i-t"^ ^ — 



□ 



Finally we give the proof of Theorem 13.71 
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Proof (Theorem I3.7p . Here we define 11 as the affine transformation such that 



u(x + n-^I^Tu + n-^l'^V + 2 ) =Tu + V + z 



for any T G R, V S 7i n and z e Ji^ . It is a consequence of Lemmas 14.11 14.31 and (j4.19p in 
Proposition 14.111 that, conditionally on Z = Tu + V G B{Q,R), the variable Z converges in law 
towards the probability measure on 7i n B{0, R) with density proportional to 

e 5 ll'^ll 3 — * (4.22) 

with respect to the Lebesgue measure on 7i n -6(0, R). Here again, the variable Z is tight thanks 
to Lemmas 14. 1[ 14.31 and Proposition 14.111 

The probability measure on tu + v ^ H, v ± u with density (|4.22p has a simple structure. It 
is clear that T and V are independent. The vector Y is determined by 

V=iV3,...,V,) 

which has a density proportional to 

"""'di;3---dwg = e-^'*^*'d^;3---d?;, 

where 

q 

Thus V is the centered Gaussian vector with covariance matrix 

= 2(<7-l)(<7-2) ('"--^^-^ 
The covariance matrix for V is computed according to Vi = and V2 = — J2i=3 ■ '-' 



(.;-l)(g-2) 11-112 



4.4 Consequences on the random-cluster model. 



Here we give the proofs of Corollary 13.81 and Proposition [ 
Proof (Corollary 13. 8p . When g > 2 is an integer, at the critical point {13, h) = (/?c,0) there are 
q + I minimizers for the free energy fjs^h, which are, on the one hand, the symmetric state 



X 



1 1 

7 ' ' ' 5 

q q 



and on the other hand, the q permutation x°-'^ of the asymmetric state 

- 1 1 1 



q q{q-'^) 9(9-1) 

We prove now that the probability of having a giant component in G{n,pn, q) has the same limit 
as the probability 

f^l3^,h„,n(— iB{x\e)] (4.23) 



for small enough e > 0, for /3„ satisfying p„ = 1 — exp(— /3„/n) and /i„ = 0. 

Indeed, let us fix a realization of the spins. Then we open edges between spins of equal color with 
probability p„ , resulting in a collection of g Erdos-Renyi random graphs G{Ni,pn, 1) for i = 1, . . . , g. 
It is known that a giant cluster appears in such a graph when lim„ NiPn > 1 (see for instance [14J). 
Yet, in the symmetric state one has lim„ NiPn = Pc/q < 1 as g > 2, hence no giant component 
appears. In the asymmetric state a;"'* on the opposite, one has limpnNi — (3c{q — l)/9 > 1 thus a 
giant component emerges with conditional probability going to 1. 

Finally, the quantity ()4.23p is computed using Theorem 13.51 after we remark that 

n \ 2 1 \n 
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□ 

Let us conclude on the computation of the partition function for the random-cluster model: 
Proof (Proposition l3 . 9p . We begin with a computation that permit to relate the partition function 
of the Curie- Weiss Potts model to that of the random-cluster model, defined at ()3.8p . Now we say 
that an edge configuration lo S {0, and a spin configurations cr G {1, . . . , 9}" are compatible 
when We = 1 = (Tj, for all e — S E{Kn), which we denote as w ^ cr. The factor q'^^") 

can be understood as the number of spin configurations cr that are compatible with lo. Hence: 

= E E n^^^^d-^)^-"^ 

ae{l,...,q}" Lje{0,l}^(^"):i^^o-eei<'„ 

E E (1-^-.,-.) 

E ^""p + ^ E '5---. 

= Z^,o,„exp(^-^(n-l)^ . (4.24) 

for (3 such that p = 1 — exp(— /3/rt). Remains to determine the asymptotics of Zp^^^n for /? < (3c- 
Thanks to the assumption (3 < [3c the minimizer of the free energy is unique and symmetric: 



X' = 



1 1 

? ' ' ' 1 

q q 



This imphes d„ = (see Remark [3^. thus TV = na;"* + ■n}/'^W. Equation (|4.13p . in the fimit 
R ^ CO, gives 

Zp„,o,n = (l-Ko„(l))A0(a^^)e-"^''"-"(^°) / e-^«-^'^('"Mw;2...d«;g (4.25) 



n 



as /p di(72 ■ • • dwq = 1 for It) = (— X]?=2 ^2, • • ■ , Wq). According to (|4.17|) one has 

Jh 

On the other hand, the free energy is easily computed: 



e-^ [1- - 



leading to 



9-1 

^/3„A« - (l + o„(l))(l-^) ' 5"exp(^^-^ 



Then (|4.24p for p„ = 1 - exp(~/3„/n) gives: 

-^p^^g.n = ^/3„,o,« exp \-^{n - 1) 



= (l + o„(l)) 1-^ <j"exp 

\ qj V 2 V 9 

and the proof is over as p„ = /3/n -f + o(l/n^) imphes = /3 -|- (7 -I- (3^/2)/n + o{l/n). □ 
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